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Abstract 

A self-consistent , non-perturbative scheme of approximation is proposed for 
arbitrary interacting quantum systems by generalization of the Hartree method. The 
scheme consists in approximating the original interaction term XHj by a suitable 
'potentiaF XV(4>) which satisfies the following two requirements: (i) the 'Hartree 
Hamiltonian' H Q generated by V{<f>) is exactly solvable i.e, the eigen states \n > and 
the eigenvalues E n are known and (ii) the 'quantum averages' of the two are equal, 
i.e. < n\Hj\n > = < n|V(0)|n > for arbitrary 'n ' . The leading-order results for 
|n > and E n , which are already accurate, can be systematically improved further 
by the development of a 'Hartree-improved perturbation theory'(HIPT) with H Q 
as the unperturbed part and the modified interaction: Afí' = X(Hj — V) as the 
perturbation. The HIPT is assured of ràpid convergence because of the 'Hartree 
condtion' : < n\H'\n >= 0. This is in contrast to the naive perturbation theory 
developed with the original interaction term XHj chosen as the perturbation, which 
diverges even for infinitesimal A ! The structure of the Hartree vacuum is shown to 
be highly non-trivial. Application of the method to the anharmonic-and double-well 
quartic-oscillators, anharmonic- sextic and octic- oscillators leads to very accurate 
results for the energy levels. In case of Xcj) 4 quantum field theory, the method 
reproduces, in the leading order, the results of Gaussian approximation, which can 
be improved further by the HIPT. We study the vacuum structure, renormalisation 
and stability of the theory in GHA. 



1. Formulation 

Consider a genèric Hamiltonian describing an arbitrary interacting quantum system: 

H = Hs + XHrtf), (1) 

where H s is exactly solvable, XHj((j>) is the interaction and À is the strength of interaction. 
The generalized Hartree approximation (GHA) to the above Hamiltonian is: 

H = H s + \V{<t>) (2) 

where V{4>) = ' Hartree potential ' (HP) [ 1 ] is the approximation to the original inter- 
action Hj((f)) and is required to satisfy the following two conditions: (i) H Q ( = ' Hartree 
Hamiltonian '(HH)) is exactly solvable, i.e., 

H Q \n > = E n \n >, < m\n >= 5 mn , (3) 

with the eigen-spectrum known and (ii) the 'quantum average'(QA) of V equals that of 
Hi,[l\ i-e, 

<n\V((j))\n> = <n\H I ((t))\n> (4) 

(The quantum average of any operator < A(cf>) > is henceforth denoted by the notation : 
< A >=< n|A|n > ). Equations (2)-(4) defining the General Hartree Approximation 
(GHA) are hence forward referred to as "Hartree Conditions (HC)". The self-consistency 
of the procedure is implicit in eqs.(2-4): the states \n > which are obtained as the solution 
of Ho, are used to specify V((f>) which, in turn, defines H . The leading approximation 
consists in finding the spectrum: \n > and E n . It may be emphasized that even the 
leading- order results capture the dominant contríbution of the interaction through the 
requirement, eq.(4), while one always deals with an exactly solvable Hamiltonian. 
Note that, because of eq.(4), the QA of H' vanishes: 

<n\H'\n>=0 (5) 

where H ' = Hj - V. 

This naturally suggests a scheme of improved perturbation theory (= 'Hartree-improved 
perturbation theory' (HIPT)) in which H is used as the unperturbed part and XH' is 
considered as the perturbation (see, below). 



2. Applications: (Quantum Mechanics) 
2. (a) The quartic - oscillator 

Consider the Hamiltonian for the (quartic) anharmonic oscillator (AHO) and the 
double-well oscillator(DWO) : 

H = \v 2 + \g<p 2 + \<P\ (6) 

where g > (< ) refers to the cases of the AHO(DWO) respectively. These 
systems have been widely studied in the literature owing to their theoretical importance 
as well as practical applications. The free-field <fi(t) and the conjugate momentum p(t) 
are parametrised in the Standard manner: 

m = ^+^i=p m = VI (6Í ~ 6) ' (7) 



/2u 

satisfying the usual commutaion relation [b,tf] = 1. Here, u, a are real constants with 
u > and a has the significance of the vacuum expectation value (VEV) of <p: o =< <fi >. 
It is to be noted that free-field operators defined by : 

1 i 
(f) = a + -j={a + a ] ), p= -^=(a ] - a) (8) 

satisfy identical commutaion relation : [a,a^] = 1. Hence (a,a^) and (b, tf) must be 
related by Quantum canonical Transformation i.e.(Boguliobov-like Transformation). De- 
noting the eigen states of the number operator b% as \n > we have the Standard definition 
\n > = (W) n \vac > j\fn\ where b\vac >= ; Wb\n >= n\n > and < m\n >= 5 rnn . 
It is to be noted that free field vacuum 10 > is defined as : a |0 0|at . The 
QA are easily calculated: < <fi >= a ; < <fi 2 >= a 2 + Ç/u; < 3 >= a 3 + 3aÇ/u and 
< 4 >= a 4 + 6a 2 Ç/u + (3/(8cj 2 ))(1 + 4^ 2 ) ; where f = n + 1/2. 

To set up the HH for the above system we propose the genèric ansatz (vàlid for 
arbitrary anharmonicity) : 

V(<j>) = A<j) 2 -B<p + C (9) 

which is amenable to immediate implementation of the HC. Self-consistency of the pro- 
cedure is built in by imposing the constraint that the constants A,B,C are chosen as 



suitable functions of < ra > . For the case of quartic anharmonicity, we determine them 
as follows: 

A =6a 2 + 3/(0/w; B = (1 + g) (tuj 2 /\ + 4 cu 2 a 3 + 12 u o f ; 
C=<0 4 > - A<(j) 2 >+B<0> (10) 

where /(£) = £ + (l/4£) . Note that, we have imposed the additional constraint of 
reproducing the variational " Gap- equation" and the "equation for the ground state" in 
the Gaussian-approximation [ 2 ], in determining the constants: A,B,C. With the choice 
above, the HH is reduced to the diagonalisable-structure, corresponding to a shifted - 
harmònic oscillator: 

Ho = ho + ^^-xf + ^p 2 , (11) 

where Q 2 = 2\A + g, x — XB/Vt 2 and h = \C — \£l 2 x 2 - However, demanding the 
consistency of the diagonalisation of H to the defining ansatz eq.(8), requires that we 
make the identification : 

to 2 = n 2 ; a 2 = x 2 , (12) 
which, in turn, imply the following constraints on u and a : 

lü 3 - u (12Xa 2 + g) - 6A/(0 = (13) 

a (AXa 2 + g+ (12\Ç/w)) = (14) 

(Here inafter,eqs.(13 ) & (14 ) are referred to as the "gap-eqn (GE)" and the "equation 
for ground state" ( EGS ) respectively). (As stated earlier, g > (< 0) refers to the case 
of the AHO(DWO)). In the case of AHO the a — 0, solution of eq(14 ) corresponds to 
the physical solution. This leads to the simplified GE, ( g = 1): 

cü 3 -cü-6À/(O=0. (15) 

This gap-equation has been derived earlier by several authors [3 - 6 ] from a variety of 
different considerations. With the aid of eqs.(ll) & (15), the energy-levels of H are given 
by 

E n = kscü + -), (16) 

4 Lü 



where 'u' is the solution of eq.(15). The numerical results for E n are given in Table — 1. 
The result is already accurate with error ~ 0.2 % to 2 % over the range of V & 'A' 
shown. 

For the case of the DWO, there are two quantum phases: ( i ) 4Àa 2 = — g — (12A£/u;), 
that leads to the spontaneously broken symmetry (SSB) phase and (ii) o — 0, that leads 
to Symmetry- Restored phase ( SR - phase ) (see,eqn.( 14 ). The former exhibits the 
double-well structure of the 'effective-potential'(for which case, the VEV of the field is 
non-vanishing: o 2 ^ 0) whereas the latter corresponds to the dynamical restoration of 
the single-well shape, with o — 0. 

The gap equations in respective phases are easily obtained : for the SSB phase ( see 
eq.(13) ) 

ul + 2^w + 6Ap(0 = (17) 
where = (5^ — 1/4^) with the physical solution given by : 

= 2yfÇ^-)c08 [\ + (^«^(A) ] (18) 

A c (0 = (-2g/3f 2 /3p(0 (19) 

where A c = "critical coupling". For the ground state ( n = 0, g= - 1, £ = 1/2) : A c (l/2) ~ 
0.09007. This shows that the SSB-phase is realized for A < A c and the SR phase is 
favoured for A > A c . The energy levels in SSB -phase are given by 

Sf ^f ( 3. + A)_|_ (20) 

For the SR phase: 

cü 3 s - gu a - 6A/(0 = (21) 
The corresponding expression for the energy-levels, are given by : 

E s n R = (Ç/4)[3u a + g/u a ] (22) 

where u s is the solution of eqn.(21) and A > A c (see,Table-2). 



2. (b) Cases of higher harmonicity 

Other cases of anharmonicity are studied in analogous manner. For the anharmonic 
sextic- oscillator, the Hamiltonian is given by 

H = \p 2 + \g<t> 2 + \<t> % (23) 

where g > (< 0) refers to the case of AHO (DWO) respectively. The Hartree- 
Hamiltonian for this case is chosen as: 

H = \? + \g 2 + A V(<f>) (24) 

with the condition, 

< H > = <H > (25) 

Identical ansatz, eq.(9), for V(<t>) is assumed: V{4>) = A (f) 2 — B(f) + C. The constants 
are calculated in the analogous manner and given as: 

A = 15 a 4 + 45^(4^ + l)/4fu; + (15/8cj 2 )(4^ 2 + 5); 

B = a [(l+^)cJ 2 /A + 6cüV 4 + 60cr 2 ^+(45/4)(4e 2 + l) ]; 

C = < 6 > —A < (f) 2 > +B < (f) > (26) 

With the above choice the Hatree-Hamiltonian, H is again reduced to the diagonalisable 
form with a shiftted field 0= — a. The structure of HH is identical to eq.(ll); together 
with eq.(12): 

H = £ + ^üü^-af + ho (27) 

where u 2 = 2\A + g, u = \Bjuj 2 and h = \C — \uj 2 u 2 . The 'gap -eqn' in this case 
is given by 

u A - uj 2 (g + 30Àcr 4 ) - 45A( ( j 2 cj/20(4^ 2 + 1) - (15À/4)(4£ 2 + 5) = (28) 
whereas, the ground-state-configuration is governed by the following equation: 

a[g + 6A((T 4 + 10Ça 2 /uj + 15(4^ 2 + 1)/8cj 2 ) ] = 0. (29) 



In the case of AHO, the a = solution of eqn. (29) corresponds to the physical 
solution. This leads to the simplified "gap-eqn.": 

c/-(^ 2 -(15A/4)(4£ 2 + 5) = (30) 

By using eqs.(27) and(30) the energy levels of H are given by 

E n = í -(2u J + ^) (31) 

where 'oo' is the solution of eq(30).The numerical results for E n are given in Table-3. 

For the case of the octic-anharmonic oscillator the Hamiltonian and the HH of the 
system are given respectively by 

H = \p 2 + \g^ + \<P 8 , (32) 

and H = ip 2 + ^0 2 + A\/(0), (33) 

with identical ansatz for V(4>): 

V(4>) = A 4> 2 — B (f) + C. The constants A,B,C are analogously determined: 

A = 28a 6 + 105<r 4 (4£ 2 + Ï)/2Çüj + (105/2cj 2 ) ( j 2 (4^ 2 + 5) + 35/i(0/2w 3 ; 

B = a [(1 + g)uj 2 /\ + 8cu 2 a 6 + 168(j 4 ^ + 105cx 2 (4£ 2 + 1) + 35£(4£ 2 + 5)/u ]; 

C = <0 8 > - A < 2 > + B < > (34) 

where MO = É 3 + (7£/2) + (9/160- 

The " gap-eqn" and eqn. for "ground- state" are as follows: 
u 5 -uj 3 (g + 56ÀO- 6 ) - 105cj 2 (Ao- 4 /0 (4£ 2 + 1) - 105cjA ( t 2 (4^ 2 + 5) - 35A/i(£) = 0, (35) 

cr [ 5 + A(8(T 6 + 168cr 4 (^/w) + 105cr 2 (4f + l)/cu 2 + 35£(4£ 2 + 5)/u; 3 ) ] (36) 

For the case of anharmonic-octic oscillator the solution o = is the physical solution. 
In this case the "gap-eqn" is simplified to 

lo 5 - gco 3 - 35A/i(0 = (37) 



which leads to the energy-levels of H 

E n = d)(5u + % (38) 
o ui 

with 'o/as the solution of eqn(37).The numerical results are given in Table-4 , compared 
with earlier results [ 7 ] . 

The generality of the above method is thus apparent from the above examples of in- 
creasing anharmonicity. In each case , the interacting system has been effectively reduced 
to an exactly solvable system while preserving the inherent non-linearity (through the 
gap-eqn.& the Hartree-condtion) of the interacting theory. The resulting energy- levels 
are uniformly accurate to within a few percent of the 'exactly' computed vàlues even in 
the zeroth-order, which is chosen to reproduce the results of the gaussian-approximation 
[ 2 ] . Further physical significance of the GHA is obtained by studying the structure and 
stability of the Hartree-vacuum as discussed below. 
3. Structure and significance of the Hartree-vacuum 

Starting from the alternative expansion of the field in terms of the 'free'-field operators, 

= a + (a + aï)/V2; p = {i/V2)(a ] - a), [a, a f ] = 1, (39) 

it follows (from eqs.(7)) that the two sets of operators must be related by a Boguliobov- 
type quantum canonical transformation : 

b = a cosh(o;) — sinh(a) 

frt — a t cos h(o;) — a sinh(a), (40) 

with the two vacua related by the transformation 

\vac > = exp[(l/2) tanh(a) (aW — aa)] |0 > ; 

a = (1/2) ln(l/u) (41) 

The set of equations: (40 ,41) imply a highly non-trivial structure of the Hartree-vacuum 
analogous to the case of the super-fluid ground state [ 8 ] and the hard-sphere-Bose gas 
[ 9 ]. In particular, the free-particle number-density in the Hartree-vacuum is non-zero 



and depends strongly on the strength of interaction: 

n =< vac\a}a\vac >= sinh 2 (a) = (1/4) [u + (1/u) - 2] (42) 

such that n ~ À 1 / 3 for À >> 1 whereas n — > for À — > (as expected). The above 
result also implies an alternative interpretation of 'u', i.e; u = (1 — uj) /(1 + u) measures 
the non-trivial structure [ 10 ] of the Hartree-vacuum. Moreover, it can be shown [ 10 ] 
that the free-fíeld ground state, |0 > is highly unstable compared to the Hartree-vacuum, 
\vac > signifying that the true vacuum is much better approximated by the latter. 

4. Perturbative Improvement 

An improved perturbation theory (HIPT) can be developed by treating the Hartree 
Hamiltonian H as the unperturbed part and XH' as the perturbation. The HIPT is, by 
construction (because of eq.(4)), guaranteed to be convergent, in contrast to the ordinary 
perturbation theory (wherein the entire interaction XHj is treated as the perturbation) 
which is divergent [ 11-13 ] even for infinitesimal À. The first- order correction in HIPT 
vanishes and the 2 nd order -correction is given by the Standard formula : 

AE™ = ]T | < n\XH'\n > \ 2 /(E n - E m ) (43) 

Only a finite number of matrix-elements contribute to eq.(43 ) corresponding to m = 
n ± 2,n ± 4 for the case of the quartic AHO (DWO) . Inclusion of this correction 
significantly improves the zeroth-order results, both in case of the AHO (Table-1) and 
the DWO (Table-2), over the full range of investigation: 0.1 < À < 100 and < n < 40. 

5. À0 4 - Quantum Field Theory 

The above formalism is easily extended to À0 4 quantum field theory described by the 
Lagrangian: C = ^(d l·l (j))(d tl (p) — ^m 2 (f) 2 — À0 4 which leads to the Hamiltonian: 

H = l [(d^/dt) 2 + (V0) 2 + m 2 2 ) ] + A0 4 (44) 

The field is Fourier-expanded in terms of the 'interacting' Fock-space operators in the 
Standard manner: 

r d 3 k 

4> = o + J — [b(k)exp(— ikx) + (k)exp(ikx) ] (45) 



where, Q k {M) = {2nf2^k 2 + M 2 = (2vr) 3 2cu fc (M), [b(k),tf(q) ] = Q k {M)5 3 (k - q), 
M = bare mass of the physical particle and \vac > is defined by b(k)\vac >= 0. 
Similarly, the physical one-particle state is defined by |k >= tf(k)\vac > . The Hartree- 
Hamiltonian is constructed in an analogous manner: 

n = \ m/dtf + (V0) 2 + m 2 2 ) ] + \V{4>) (46) 

with V{4>) parameterised as before : 

V(<j>) = A<f) 2 - B<p + C (47) 

The various QA's are now given by : < 4> 2 >= o 2 + I , < 4 >= ct 4 + 6a 2 I + 3/q, 
< (d(j)/dt) 2 >= h, < (V0) 2 >= h - M 2 I , where < À >= < vac\À\vac > and 

/„ = /„(M 2 ) = I [col(M) f, n = 0, ±1, ±2, ±3, (48) 

are the Stevenson-Integrals [ 14 ]. Applying the Hartree-condition , eq.(4), and ensuring 
equivalence to the GEP [ 14 ] leads, as before, to the complete determinat ion of the 
parameters in V(4>): 

A = 6 < 2 > ; B = 8 a 3 and C = 3a 4 - Qa 2 I - 3/ 2 (49) 

By analogous procedure, the Hartree - Hamiltonian is then reduced to diagonal form 
(corresponding to an effectively free Klein-Gordon theory) by a shift of the field and 
energy : 

Ho = \ [(d V /dt) 2 + (Vv?) 2 + MV) ] + h (50) 

where <p = - a and h = - (l/2)M 2 a 2 + A C, with 'C given in eq.( 49 ) 
The 'gap-equation' is given by 

M 2 = m 2 + 12Xa 2 + 12A/ (M 2 ) (51) 

and the VEV now satisfies the consistency condition: 

a[M 2 - 8 A a 2 } = (52) 



X R = (l/4\)d 4 U/daX=o = A 



(55) 



The physical solution of this equation for the vacuum-configuration is at o = 0, 
as intuitively expected for the symmetric À0 4 theory (i.e. with m 2 > 0) considered 
here. This is further verified by the computation of the effective potential and the 
renormalized parameters as shown below. 

Employing the Standard definition of the effective potential U(a): < vac\H\vac > = 
< vac\TLo\vac >= U(<r), we have, 

U(a) = h- 3A/ 2 + (l/2)mV + Xa 4 (53) 

It is important to note that, eqs.( 51 - 53 ) coincide with the results based on the 
Gaussian Effective Potential (GEP) [ 14 ] . It is then straight- forward to carry out 
the renormalisation program as has been done in ref.[14 ]. The renormalized parameters 
are given by the following expressions: 

m\ = d 2 U/da\ =0 = m 2 + 12À/ (M 2 ); M 2 = M 2 [o = 0) (54) 

1 - 12XI-! 
1 + 6A/_i 

It has been shown [14] in case of the GEP that a non— trivial version of À0 4 theory emerges 
which can not be realized in lattice theories (for any fínite lattice-spacing). Identical 
conclusion holds for the GHA in the leading order, because of the equivalence established 
here. 

Beyond proving the equivalence of the GHA to GEP in the leading order, the new 
results derived are : 

(i) the vacuum — structure: 

\vac >= ex P {\ j ^y^(k) [^(k^-k) - a(k)a(-k) ]}|0 > (56) 

where /3(k) is analogous to a(k) (see,eqn.(40) with a — > a(k)). The "vacuum- structure- 
function" is denoted by u(k) given by 

u(k) = cü k (m)/üü k (M) = K^ + m^/í^ + M 2 )] 1 / 2 (57) 

The free-particle number - density in the Hartree - vacuum is given by (V = spatial 
volume ) 

n(k) =< vac\a ] '(k)a(k)\vac > /VVt k (m) (58) 



One then calculates p(k) as defined below 

lim n(k)/ra(0) = p(k) = [1 + (k 2 /m|) Y 1 ' 2 (59) 

A^oo 

where n(0) = (m/m R )/(32Tc 3 ) is the maximum value of n(k). Eqs.( 57-58) show the 
condensate- structure of the physical vacuum consisting of correlated off- shell particle 
pairs ( with momenta k and -k) with a non-trivial dependence on |k|! This condensate 
- structure is expected to play significant role in the thermodynamic- properties of the 
system, if the behaviour persists to non-zero value of temprature, T. 

The static- potential U(r) of the À0 4 theory is easily calculable from the 2-particle 
correlation function C/(x — y) defined as: 

C/(x — y) =< t>ac|0(x, 0)<f>(y, 0) \vac > / < vac\vac > (60) 

which leads to 

U(r) = - j —exp^k.r)(k 2 + m 2 R )- 1/2 (61) 

where r = x — y. The above integral can be related to the modified Bessel-function 
K\{m R r) and one obtains: 

U(r) = runK^mR^/A^r, (62) 

where, r = |r|. Note that U(r) diverges in the limit r — > 0, as it should (in any local 
quantum field-theory) and behaves asymptotically as 

lim U(r) ~ r~ 3 / 2 exp(—m R r) (63) 

which shows a fall — off faster than the Yukawa — potentiah To our knowledge, the 
important information (contained in Eqs.( 60 - 63 )) regarding the static-potential of the 
symmetric À0 4 theory are new results. 
6. Summary and Conclusion 

The generalised Hartree approximation is, in principle, applicable to arbitrary quantum- 
systems with interaction. The method consists in mapping the interacting system to an 
exactly solvable system while preserving the essential non-linearity and the major effects 



of the interaction of the original system through the self-consistent feed-back mechanism. 
The bàsic approximation, which is non — perturbative in character, lends itself to be 
systematically improved by modified perturbation theory developed about the Hartree- 
vacuum and the Hartree- Hamiltonian chosen as the unperturbed part and the modi- 
fied interaction term XH' taken as the perturbation. This modified perturbation the- 
ory can be shown to be rapidly convergent in contrast to the divergence of the naive 
perturbation theory developed about the free-field vacuum with the entire interaction 
term XHi chosen as perturbation. We conjecture that the instahility of the free — 
fíeld vacuum and the divergence of perturbation theory developed about this vacuum 
may be intimately related . The structure of the Hartree-vacuum emerges to be highly 
non-trivial with a strong dependence on À analogous to the case of the ground states of 
superfluid Helium, hard- sphere Bose-gas etc. 

The method and its perturbative improvement applied to the case of the quartic, 
sextic and the octic anharmonic - oscillator and the quartic double-well oscillator leads to 
excel·lent results over the entire allowed range of the interaction strength À and excitation 
energy level 'n' . The generalisation to the case of À0 4 field theory reproduces, in the 
leading order, the results derived from the Gaussian effective potential [ 14 ]. Going 
beyond the leading order through the HIPT, the results of the Gaussian approximation 
can be systematically improved. The structure of the Hartree-vacuum again emerges quite 
non-trivial and interesting, being characterized by the condensation of off-shell, correlated 
particle-pairs and leading to the definition of a 'vacuum — structure function'. It has also 
been established elsewhere [ 10 ] that the effective potential based upon the free-field 
('perturbative') vacuum leads either to instability or to triviality. We also derive the 
inter — particle (static) potential by evaluating the two-particle correlation-function in 
the Hartree vacuum and show that the potential falls faster than the Yukawa potential! 
The generalization to finite temperature appears straight-forward. 
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Table- 1: Sample results for the zeroth - order (GHA) compared with results of earlier 
calculations from ref.[ 3 ] (shown in parentheses), over a wide range of 'A' and 'n'. Also 
shown are the results (in square brackets) of the Hartree -improved perturbation theory 
(HIPT) up to second order (see text). 



A 


Eo 


E 1 


E2 


E 4 


E10 


E40 


0.1 


0.56031 


1.7734 


3.1382 


6.2052 


17.2267 


94.84 




(0.55915) 


(1.7695) 


(3.1386) 


(6.2203) 


(17.352) 


(90.56) 




[0.55911] 


[1.7694] 


[3.1391] 


[6.2239] 


[17.374] 


[95.766] 


1.0 


0.81250 


2.7599 


5.1724 


10.902 


32.663 


192.79 




(0.80377) 


(2.7379) 


(5.1792) 


(10.902) 


(32.963) 


(194.60) 




[0.80321] 


[2.7367] 


[5.1824] 


[10.982] 


[33.013] 


[195.15] 


10.0 


1.5313 


5.3821 


10.3240 


22.248 


68.177 


409.89 




(1.5050) 


(5.3216) 


(10.3471) 


(22.409) 


(68.804) 


(413.94) 




[1.5030] 


[5.3177] 


[10.356] 


[22.457] 


[68.996] 


[415.18] 


100.0 


3.1924 


11.325 


21.853 


47.349 


145.843 


880.55 




(3.1314) 


(11.187) 


(21.907) 


(47.707) 


(147.231) 


(889.32) 




[3.1266] 


[11.178] 


[21.927] 


[47.817] 


[147.652] 


[892.03] 


1000.0 


6.8280 


24.272 


46.902 


101.742 


313.720 


1895.90 




(6.6942) 


(23.972) 


(47.017) 


(102.516) 


(- -) 


(- -) 




[6.6836] 


[23.952] 


[47.062] 


[102.75] 


[317.65] 


[1920.70] 



Table- 2: The computed energy levels of the quartic - DWO in the zeroth order of GHA 
compared with the earlier calculations including twenty - orders of perturbation theory 
[3] shown for sample vàlues of 'A' and 'n' . Also shown are the results after inclusion of 
second order correction in HIPT, denoted as E^f 1 . 



X 


n 




E^ 


Ref[3] 


0.1 





0.5496 


0.4606 


0.4702 




1 


0.8430 


0.7553 


0.7703 




2 


1.5636 


1.6547 


1.6300 




4 


3.5805 


3.7232 


3.6802 




10 


12.192 


12.517 


12.400 


1.0 





0.5989 


0.5752 


0.5800 




i 

X 


9 1 9^0 


9 nsnn 






2 


4.2324 


4.2600 


4.2500 




4 


9.4680 


9.5950 


9.5600 




10 


30.530 


30.650 


30.420 


10.0 





1.4098 


1.3752 


1.3800 




1 


5.0650 


4.9910 


5.0900 




2 


9.8660 


9.9050 


9.8900 




4 


21.561 


21.791 


21.700 




10 


66.950 


67.820 


67.620 


100.0 





3.1340 


3.0650 


3.0700 




1 


11.175 


11.024 


11.002 




2 


21.638 


21.715 


21.700 




4 


47.023 


47.505 


47.200 




10 


145.27 


147.10 


146.70 



Table- 3: Sample results for the zeroth - order (GHA) for the sextic - AHO compared 
with results of earlier calculations from ref.[ 7 ] (shown in parentheses), over a wide range 
of 'À' and 'n' . Percentage of error is shown in square bracket. 



p = 


0.2 


2.0 


10.0 


100.0 


400.0 


2000.0 


Eq 


1.193 


1.676 


2.323 


3.947 


5.521 


8.206 




(1.174) 


(1.610) 


(2.206) 


(3.717) 


(5.188) 


(7.702) 




[1.611] 


[4.079] 


[5.313] 


[6.188] 


[6.415] 


[6.544] 


Ex 


3.966 


5.931 


8.420 


14.52 


20.39 


30.37 




(3.901) 


(5.749) 


(8.115) 


(13.95) 


(19.56) 


(29.12) 




[1.681] 


[3.165] 


[3.762] 


[4.148] 


[4.244] 


[4.298] 


E 2 


7.420 


11.61 


16.74 


29.16 


41.03 


61.18 




(7.382) 


(11.54) 


(16.64) 


(28.98) 


(40.78) 


(60.81) 




[0.523] 


[0.612] 


[0.6179] 


[0.6157] 


[0.6145] 


[0.6138] 


E 4 


16.15 


26.48 


38.73 


68.01 


95.90 


143.2 




(16.30) 


(26.83) 


(39.29) 


(69.05) 


(97.38) 


(145.4) 




[0.9170] 


[1.302] 


[1.426] 


[1.499] 


[1.517] 


[1.527] 




26.88 


45.08 


66.36 


117.0 


165.1 


246.5 




(27.29) 


(45.94) 


(67.70) 


(119.4) 


(168.5) 


(251.7) 




[1.50] 


[1.870] 


[1.98] 


[2.043] 


[2.058] 


[2.067] 




53.24 


91.17 


135.0 


238.7 


337.1 


503.8 




(54.31) 


(93.26) 


(138.2) 


(244.5) 


(345.3) 


(516.1) 




[1.967] 


[2.245] 


[2.323] 


[2.367] 


[2.377] 


[2.383] 


E\A 


85.01 


147.0 


218.3 


386.6 


546.2 


816.3 




(86.78) 


(150.4) 


(223.4) 


(395.7) 


(559.1) 


(835.6) 




[2.047] 


[2.230] 


[2.279] 


[2.306] 


[2.313] 


[2.316] 


En 


111.9 


194.4 


289.0 


512.1 


723.7 


1082.0 




(114.0) 


(198.3) 


(294.9) 


(522.7) 


(738.6) 


(1104.0) 




[1.868] 


[1.974] 


[2.001] 


[2.016] 


[2.020] 


[2.022] 



Table- 4: Sample results for the octic- AHO in the zeroth - order (GHA) compared 
with results of earlier calculations from ref . [ 7 ] (shown in parentheses) , over a wide range 
of 'A' and V. 



A = 


0.1 


1.0 


5.0 


50.0 


200.0 


Eq 


1.3005 


1.7794 


2.3290 


3.5565 


4.6425 




(1.2410) 


(1.6413) 


(2.1145) 


(3.1886) 


(4.1461) 




4.4717 


6.3946 


8.5167 


13.172 


17.259 




(4.2754) 


(5.9996) 


(7.9296) 


(12.1950) 


(15.9519) 


E 2 


8.6264 


12.717 


17.126 


26.698 


35.062 




(8.4530) 


(12.421) 


(16.711) 


(26.033) 


(34.183) 


E 4 


19.763 


30.026 


40.863 


64.165 


84.444 




(19.9930) 


(30.4605) 


(41.4947) 


(65.20180) 


(85.8251) 




34.217 


52.669 


72.044 


113.48 


149.47 




(35.0560) 


(54.1403) 


(74.0830) 


(116.7629) 


(153.8278) 


Eg 


51.570 


80.013 


109.65 


172.99 


227.97 




(53.145590) 


(82.6496) 


(113.3486) 


(178.9215) 


(235.8193) 


Eg 


61.239 


95.255 


130.64 


206.23 


271.81 




(63.2253) 


(98.5529) 


(135.2598) 


(213.6157) 


(281.5864) 




71.532 


111.49 


153.01 


241.64 


318.52 




(73.9545) 


(115.4899) 


(158.5991) 


(250.5751) 


(330.3433) 


En 


824.24 


128.68 


176.69 


279.14 


368.00 




(85.3079) 


(133.4201) 


(183.3103) 


(289.7106) 


(381.9720) 


E\i 


93.893 


146.79 


201.65 


318.67 


420.14 




(97.2636) 


(152.3080) 


(209.3443) 


(330.9440) 


(436.3695) 


E13 


105.92 


165.79 


227.84 


360.14 


474.85 




(109.7967) 


(172.1125) 


(236.6436) 


(374.1834) 


(493.4143) 


Eia 


118.49 


185.65 


255.21 


403.50 


532.06 




(122.8909) 


(192.8082) 


(265.1732) 


(419.3737) 


(553.0335) 



